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Periconceptional folic acid use can often prevent neural tube defects (NTDs). Variants of 
genes involved in folate metabolism in mothers and children have been associated with 
occurrence of NTDs. We identified Irish families with individuals affected by neural tube 
defects. In these families, we observed that neural tube defects and birth defects overall 
occurred at a higher rate in the maternal lineage compared with the paternal lineage. 
The goal of this study was to look for evidence for genetic effects that could explain 
the discrepancy in the occurrence of these birth defects in the maternal vs. paternal 
lineage. We genotyped blood samples from 322 individuals from NTD-affected Irish 
families, identified through their membership in spina bifida associations. We looked for 
differences in distribution in maternal vs. paternal lineages of five genetic polymorphisms: 
the DHFR 19 bp deletion, MTHFD1 1958G>A, MTHFR 1298A>C, MTHFR 677C>T and 
SLC19A1 80A>G. In addition to looking at genotypes individually, we determined the 
number of genotypes associated with decreased folate metabolism in each relative 
("risk genotypes") and compared the distribution of these genotypes in maternal vs. 
paternal relatives. Overall, maternal relatives had a higher number of genotypes associated 
with lower folate metabolism than paternal relatives (p = 0.017). We expected that 
relatives would share the same risk genotype as the individuals with NTDs and/or their 
mothers. However, we observed that maternal relatives had an over-abundance of any 
risk genotype, rather than one specific genotype. The observed genetic effects suggest 
an epigenetic mechanism in which decreased folate metabolism results in epigenetic 
alterations related to the increased rate of NTDs and other birth defects seen in the 
maternal lineage. Future studies on the etiology of NTDs and other birth defects could 
benefit from including multigenerational extended families, in order to explore potential 
epigenetic mechanisms. 

Keywords: neural tube defects, folate metabolism, DHFR 19bp deletion, MTHFD1 1958G>A, MTHFR 1298A>C, 
MTHFR 677C>T, SLC19A1 80A>G, maternal inheritance 



INTRODUCTION 

Neural tube defects (NTDs) are congenital anomalies of the brain 
and spinal cord that cause significant morbidity and mortality 
(Sutton et al., 2008). Historically, Ireland has had high rates of 
NTDs (Coffey, 1983), with rates decreasing significantly in the 
last 50 years (Botto et al., 2006). Both environmental and genetic 
factors have been implicated in the etiology of NTDs. The evi- 
dence for environmental factors lies in the protective effect of 
periconceptional folic acid intake (Czeizel and Dudas, 1992). The 
evidence for genetic factors lies in the increased recurrence risk in 
siblings of individuals with NTDs, the increased risk for NTDs 
among other relatives, and the increased risk for birth defects 
overall among both close and distant relatives (MRC Vitamin 
Study Research Group, 1991; Deak et al., 2008). Maternal relatives 



have higher NTD and overall birth defect rates than paternal rel- 
atives, consistent with a combination of genetic and epigenetic 
mechanisms contributing to NTD and overall birth defects risk 
(Byrne et al, 1996; Byrne, 2008, 2010, 201 1; Deak et al, 2008). 

The genetic mechanisms associated with NTDs have been 
addressed through studies on associations between NTD risk and 
variants in genes involved in folate metabolism, chosen because 
of the protective effect of folic acid on NTDs (Molloy et al., 
2009; Shaw et al., 2009). Methylene tetrahydrofolate reductase 
(MTHFR) plays an important role in folate metabolism. The 677 
T allele results in reduced enzyme activity and decreased folate 
metabolism (Frosst et al, 1995). The MTHFR 677 TT geno- 
type in infants and mothers has been associated with increased 
risk for NTDs [reviewed in Botto and Yang (2000)]. Studies 
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in Irish populations found significantly more individuals with 
NTDs with the TT and CT genotypes and more parents of 
affected individuals with the TT genotype (Whitehead et al., 1995; 
Shields et al., 1999; Kirke et al, 2004; Carter et al, 2011). The 
MTHFR 1298A>C variant also results in reduced enzyme activity, 
although its role in NTD risk is less clear (van der Put et al., 1998; 
Botto and Yang, 2000; Wang et al, 2012a). Parle-McDermott et al. 
(2003) did not find an association between the MTHFR 1298A>C 
variant and NTDs in an Irish population. Methylene tetrahydro- 
folate dehydrogenase 1 (MTHFD1) acts in folate metabolism and 
DNA synthesis, and the 1958G>A variant results in decreased 
enzyme stability (Christensen et al., 2009). Studies in Irish pop- 
ulations (Brody et al, 2002; Parle-McDermott et al, 2006) iden- 
tified a significant excess of MTHFD1 AA homozygotes among 
mothers of NTD cases compared to controls. A recent larger scale 
study in an Irish population also identified the MTHFD1 A allele 
as a risk factor for NTDs (Pangilinan et al., 2012). Dihydrofolate 
reductase (DHFR) also acts in folate metabolism. Studies have 
yielded conflicting results on the effect of the 19 bp deletion 
at intron 1 on red blood cell folate concentrations (Kalmbach 
et al, 2008; Stanislawska-Sachadyn et al., 2008). Researchers from 
Trinity College in Dublin observed a slight protective effect of 
the homozygous DHFR deletion genotype in mothers (Parle- 
McDermott et al., 2007). SLC19A1 plays an important role in 
folate transport, and some previous studies have reported associ- 
ations between the 80A>G variant and NTD risk (Chango et al., 
2000; De Marco et al, 2003; Wang et al, 2012b). O'Leary et al. 
(2006) did not find a difference in the frequency of the SLC19A1 
80A>G reduced folate carrier in mothers, fathers and affected 
individuals compared with controls in an Irish population. 

Epigenetic mechanisms related to NTDs have not been well 
explored. The higher rates of NTDs in maternal relatives com- 
pared with paternal relatives suggest the presence of epigenetic 
effects. Folate metabolism genes are of interest because folate 
metabolism is critical to the production of one-carbon sub- 
units needed for epigenetic modification (Crider et al, 2012). 
Decreased folate metabolism, seen with certain gene variants, 
might hinder proper establishment or maintenance of epigenetic 
modifications. The epigenetic changes of interest would occur in 
the maternal genome in utero, that is, during the grandmother's 
pregnancy (Crider et al., 2012). Thus, in contrast to genetic mech- 
anisms, in which the mothers' and affected individuals' genotypes 
would be most relevant, for epigenetic mechanisms, the mater- 
nal and paternal grandparents' genotypes would be important. 
In the absence of information on grandparents' genotypes, aunts, 
uncles, and first cousins provide an indirect measure of the 
grandparents' genotypes. 

In this hypothesis-generating study, we investigated possible 
mechanisms for the excess risk for NTDs and birth defects over- 
all in the maternal lineage of Irish individuals with NTDs. To do 
this, we evaluated five polymorphisms in genes involved in folate 
metabolism in the maternal vs. paternal lineage, both individually 
and with low folate-metabolizing genotypes grouped together. To 
our knowledge, this study is the first to look for genomic asso- 
ciations in extended relatives in NTD families and the effect of 
maternal vs. paternal lineage risk. We chose five polymorphic 
markers for evaluation, which were among the most intensely 



studied so far: the 19 bp deletion in DHFR, MTHFD1 1958G>A, 
MTHFR 1298A>C, MTHFR 677C>T, and SLC19A1 80A>G. 

MATERIALS AND METHODS 
STUDY DESIGN AND PARTICIPANTS 

The participants in this cross-sectional study were members 
of Irish families with NTDs. These families had previously 
participated in epidemiologic studies of the Boyne Research 
Institute [for instance, (Byrne, 2011)]. The families were iden- 
tified initially through their membership in the Louth-Meath 
(northeastern region of the Republic of Ireland) branch of the 
Irish Association for Spina Bifida and Hydrocephalus, and the 
Ballymena (Northern Ireland) branch of the (UK) Association 
for Spina Bifida and Hydrocephalus. Families were also recruited 
through word of mouth and through radio announcements in the 
Louth-Meath area. 

To be eligible for the study, male and female participants had 
to be 15 years of age or older, reside in the northeast region of 
the Republic of Ireland or Northern Ireland, and be either the 
individual with an NTD (only one NTD in sibship) or a blood rel- 
ative of the affected individual. In terms of relationship, relatives 
(defined according to their relationship to the affected individual) 
had to be one of the following: parents, siblings, nieces, nephews, 
offspring (of the affected individual), grandparents, uncles, aunts, 
first cousins and first cousins once removed (the offspring of the 
first cousins). 

Family members were interviewed in phases. Phase I evalu- 
ated pregnancy histories and the health of siblings in the original 
nuclear families (1995-2002). Phase II subjects (uncles and aunts 
of the affected individuals) were identified by their siblings during 
Phase I and interviewed between 2000 and 2002. Phase III con- 
sisted of interviews with first cousins of the affected individuals 
(the children of the uncles and aunts) between 2002 and 2009. 
Between 2005 and 2006, affected individuals and their siblings 
were interviewed. 

For the present study, conducted in 2007, the response rate 
among eligible, locatable relatives was 49.5%, or 331 participants. 
From these, 322 gave a blood sample. Of these with a blood sam- 
ple, 34 were relatives in families in which the individual with 
an NTD had anencephaly, 41 in which the affected individual 
had spina bifida occulta, 240 in which the affected individual 
had open spina bifida and 7 in which the affected individual had 
encephalocele. 

This study was conducted according to the guidelines from 
the Declaration of Helsinki and all procedures involving human 
subjects were approved by the Ethics Board of the Boyne 
Research Institute and the Ethics Committee of the Health Service 
Executive of the Irish Government's Department of Health 
and Children. Written informed consent was obtained from all 
subjects. 

STATISTICAL ANALYSIS 

Data were managed using Microsoft EXCEL®, and statistical anal- 
ysis was conducted using SAS 9.1 (SAS Institute, Cary, NC, USA), 
STATA®, and IBM® SPSS® vl9. Measures of association were 
evaluated with f-tests for continuous variables and chi-square 
and Fisher's exact tests for categorical variables, with the level of 
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statistical significance set at 0.05. First-degree relatives were par- 
ents, siblings, and children; second-degree relatives were aunts, 
uncle, nieces, nephews, and grandparents; and third-degree and 
higher relatives were cousins and first cousins once removed. 
Genotypes were analyzed using an additive model, consistent 
with previous studies. For logistic regression analyses, genotypes 
were dichotomized using a recessive model, consistent with pre- 
vious studies. For each individual, we calculated the number of 
"risk genotypes" present. "Risk genotypes" were those associated 
with lower folate metabolism in previous studies (DHFR WT/WT, 
MTHFD1 1958 AA, MTHFR 1298 CC, MTHFR 677 TT, and 
SLC19A1 80 GG). 

DNA METHODS 

DNA was extracted from blood clots using Centra Puregene 
Blood kit and Clotspin Baskets per manufacturer's protocol 
(Qiagen, Valencia, CA). Quantitation was performed by spec- 
trophotometric reading at 260nm and all samples normalized 
to 10ng/u,l. For PCR reactions, 2 u,l genomic DNA (20 ng) 
was used. Pyrosequencing technology (Qiagen, Valencia, CA) 
was used to genotype DNA samples for DHFR 19bp dele- 
tion (rs70991108), MTHFD1 1958G>A (rs2236225), MTHFR 
1298A>C (rsl801131), MTHFR 677C>T (rsl801133), and 
SLC19A1 80A>G (rsl051266). PCR and pyrosequencing primers 
(Table 1 ) were designed using PSQ Assay Design Software Version 
1.0.6. After amplification, the pyrosequencing primer was used 
to extend the PCR product by addition of nucleotides in a spe- 
cific sequential order based on the known sequence. The resulting 
pyrogram of signal peaks was used to determine the genotype. 
PCR amplification for the DHFR 19bp deletion was performed 
using the GC Rich PCR System (Roche, Indianapolis, IN) as 
previously described (Parle-McDermott et al., 2007). PCR for 
all other variants was performed in a 20 \il reaction mix con- 
taining 0.5 u,M each forward and reverse primer and Applied 
Biosystems GeneAmp® Fast PCR Master Mix. Amplification was 
performed using the 9800 Fast PCR System according to man- 
ufacturer's suggested protocol (Applied Biosystems, Foster City, 
CA). Pyrosequencing was carried out as described by manufac- 
turer's protocol (Qiagen, Valencia, CA). 

QUALITY CONTROL 

Five percent of the samples were selected at random for repeat 
testing of the five polymorphisms, and the concordance was 
checked. We also tested each variant for any significant departure 



from Hardy- Weinberg equilibrium (HWE) by a chi-square 
goodness-of-fit test to verify data quality. In addition, Mendelian 
inheritance for all variants was assessed in 13 families with full 
trios using the PedCheck program (O'Connell and Weeks, 1998). 
Water and genotype controls (Coriell Cell Repositories, Camden, 
NJ) were included on each assay plate. The tests for HWE indi- 
cated that none of the variants deviated significantly from HWE, 
supporting high quality genotype data. All replica test results were 
100% concordant. 

RESULTS 

We obtained blood samples from all categories of relatives, 
including individuals with NTDs (n = 13), offspring (N = 2), 
fathers (N = 29), mothers (N = 37), siblings (N = 59), nieces 
and nephews (N = 2), uncles and aunts (N = 66), first cousins 
(N = 105), and grandparents (N = 2) from 47 families. Paternal 
relatives (N = 65) and maternal relatives (N = 115) comprised 
uncles and aunts, first cousins, first cousins once removed, and 
grandparents, who were related to the affected individual through 
the father and mother, respectively. The average age of the study 
participants was 41 years. The overall distribution of each of the 
five polymorphisms is shown in Table 2. The genotype distribu- 
tions for the five polymorphisms in our study were similar to 
those from other studies, with a few exceptions (Table 3). 

ANALYSIS OF INDIVIDUAL POLYMORPHISMS 

The distributions of each genotype and allele by relationship to 
the individual with an NTD, degree of relationship to the affected 
individual, and maternal vs. paternal lineage are shown in Table 4. 
We did not observe any significant associations between individ- 
ual genotypes and the affected individuals or mothers, although 
the higher folate-metabolizing MTHFR 677 CC genotype tended 
to be found less often in affected individuals (15.4%) than in all 
other relatives (42.7%) (Fisher's exact test, p = 0.08) (Table 4 and 
data not shown). We also did not see any significant associations 
between individual genotypes and maternal vs. paternal lineage 
(Table 4). 

Relatives of individuals with NTDs might be expected to 
have genotypes associated with increased or decreased risk for 
NTDs, and the prevalence of these genotypes could increase or 
decrease with degree of closeness to the affected individual. For 
MTHFR 1298A>C, a lower percentage of first-degree relatives 
(34.6%) had the higher folate-metabolizing AA genotype com- 
pared with second-degree (50.0%) relatives and third-degree and 



Table 1 | Primer sequences and PCR amplicon size for DHFR, MTHFD1, MTHFR, and SLC19A 7 variant analysis. 



Variant 


Forward primer (5'-3') a 


Reverse primer (5'-3') a 


Sequencing primer (5'-3') b 


Amplicon 
size (bp) 


DHFR 19bp deletion (rs70991108) 


CAAGAACGGGGACCTGCC 


•ATCCTCTCGCCGGGAGTC 


GCCCACGGTCGGGGT 


114,133 


MTHFD1 1958G>A (rs2236225) 


TCATTTTG CTTAG AG CACAGTAG A 


* CATC C ACTTTC CTG G CTTAAAT 


CACAGTAGAGAGTGCCAA 


63 


MTHFR 1298A>C (rs1801131) 


G G AG CTG CTG AAG ATGTG G 


•TGGTTCTCCCGAGAGGTAAAG 


AGGAGCTGACCAGTGA 


85 


MTHFR 677C>T (rs1801133) 


•ATTGGCAGGTTACCCCAAAGG 


TCGGTGCATGCCTTCACAA 


CGTGATGATGAAATCG 


156 


SLC19A1 80A>G (rs1051266) 


•GTGGAACCTGGGCCTGAC 


GTAGGGGGTGATGAAGCTCTC 


CAAAG GTAG CACACGA 


114 



3 'Biotinylated at 5 end. 

b Represent forward primers for DHFR, MTHFD1, MTHFR 1298, and reverse primers for MTHFR 677 and SLC19A1. 
c MTHFR 677 and SLC19A 7 assays on reverse complement strand. 
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Table 2 | Characteristics of the participants, northeast region of the 
Republic of Ireland and Northern Ireland, 1995-2009. 







Range 


Median 


Age in years 


Overall 


16-86 


41.5 




Affected individual 
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N 
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Gender 


Males 


120 
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Females 


202 
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Affected individual 
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Z 


u.o 
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zy 


y .U 




iviomer 
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Sibling 


oy 
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z 
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Uncle/aunt 


66 


20.5 




First cousin 


105 


32.6 




First cousin once removed 


7 


2.2 




Grandfather 


1 


0.3 




Grandmother 


1 


0.3 




1st degree 


127 


39.4 




2nd degree 


70 


21.7 




3rd degree or higher 


112 


34.8 


Genotyping results 


Unrn iyop aeietion 








Hal/Hal 


O I 


10 n 
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1 O.O 
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A A 

AA 


145 


45.0 




AL- 


14/ 


A C "7 
4o. / 






on 
oU 


Q Q 

y .0 




MTHFR 677C>T 








CC 


134 


41.6 




CT 


156 


48.5 




jj 


32 


9 9 




SLC19A1 80A>G 








AA 


58 


18.0 




AG 


168 


52.2 




GG 


96 


29.8 


Line 


Paternal relatives 


65 


36.0 




Maternal relatives 


115 


64.0 



higher (53.6%) relatives (p = 0.009). Similarly, the A allele was 
less common in first-degree relatives (61.4%) than in second- 
degree relatives (70.0%) and third-degree and higher (73.7%) 
relatives (p = 0.005). When the genotype (AA vs. A/C, C/C) 



and relatives (first-degree relatives vs. all other relatives) were 
dichotomized, the odds ratio for the association of the AA geno- 
type with first-degree relatives was 0.49 (0.30-0.78, p = 0.002) 
(data not shown). Consistent with this, a lower percentage of sib- 
lings, 30.5%, had the AA genotype compared to 48.3% among 
other relatives (p = 0.04, Table 4). We did not observe a consis- 
tent trend of genotypes being more or less common in relatives 
more closely related to the affected individual for any of the other 
genotypes (Table 4). 

We observed some associations with individual categories of 
relatives. For example, fewer first cousins had the higher folate- 
metabolizing DHFR del/del genotype compared with all other 
relatives (12.4 vs. 22.1%, p = 0.11, Table 4 and data not shown), 
the higher folate-metabolizing MTHFD1 GG genotype was higher 
in uncles and aunts compared with all other relatives (51.5 vs. 
31.6%, p = 0.003) (Table 4), and the lower folate-metabolizing 
SLC19A1 GG genotype was more common in uncles and aunts 
than in other relatives (42.4 vs. 26.6%, p = 0.04) (Table 4). 
However, the meaning of these isolated associations is unclear. 

NUMBER OF RISK GENOTYPES 

While the frequencies of individual genotypes did not differ sig- 
nificantly between maternal and paternal relatives (Table 4), the 
number of risk genotypes, defined as the number of homozy- 
gous genotypes associated with lower folate metabolism, did 
show significant variation between the maternal and paternal 
lines (p = 0.02) (Table 5). A lower percentage of maternal rela- 
tives (20.9%) had no (0) risk genotypes compared with paternal 
relatives (41.5%), while more maternal relatives (54.8%) than 
paternal relatives (35.4%) had 1 risk genotype. When the risk 
genotypes were dichotomized (0, 1-3), the odds ratio for associ- 
ation with the maternal line compared with the paternal line was 
2.69 (1.38-5.25,/?= 0.004). 

The numbers of risk genotypes in affected individuals and 
mothers (p = 0.004) or mothers alone (p = 0.005), compared 
with all other relatives were also statistically significant (Table 5). 
A higher percentage of affected individuals and mothers (30.0 
vs. 17.3%) and mothers alone (32.4 vs. 17.5%) had 2 risk geno- 
types compared with all other relatives, while a lower percentage 
of affected individuals and mothers (22.0 vs. 47.8%) and mothers 
alone (18.9 vs. 47.1%) had 1 risk genotype. A higher percentage 
of affected individuals and mothers (44.0 vs. 31.6%) and mothers 
alone (46.0 vs. 31.9%) had 0 risk genotypes. 

DISCUSSION 

In the Irish families in our study, we found that maternal relatives 
were more likely than paternal relatives to have one or more geno- 
types associated with decreased folate metabolism. As expected 
from previous studies, mothers and affected individuals had dif- 
ferential distributions of risk genotypes compared with all other 
relatives. 

The observed increased number of affected relatives in the 
maternal lineage of families could be due to genetic and/or epi- 
genetic mechanisms. A simple genetic explanation would be that 
more maternal relatives have alleles associated with increased risk 
for NTDs simply due to chance, so the likelihood of inheriting 
these alleles from the maternal side would be increased. With a 
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Table 3 | Distribution of genotypes and comparison with other studies. 
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Parle-McDermott et al., 


2007 


Ireland, controls 


46 




18.6 


128 




51.8 


73 




29.6 
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Kalmbach etal., 2008 




US, unaffected 


216 




17.1 


646 




51.4 


396 




31.5 


0.39 


This study 




Ireland, mothers 


7 




18.9 


22 




59.5 


8 




21.6 
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2007 


Ireland, mothers 
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48.6 
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54 




47.8 


39 




34.5 


0.90 
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2006 


Ireland, controls 
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19.5 
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53.4 
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27.1 


0.02 


This study 




Ireland, Mothers 


8 




21.6 


17 




46.0 


12 
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Parle-McDermott et al., 
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65 
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62 
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30 
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Utah 


49 




43.4 
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Ireland, controls 


519 




49.4 
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O'Leary et al., 2006 




Ireland, controls 


103 




27.0 


179 




47.0 


99 




26.0 


0.02 


This study 




Ireland, affected individuals 


4 




30.8 


7 




53.9 


2 




15.3 


Ref. 


O'Leary et al., 2006 




Ireland, affected individuals 


61 




23.0 


131 




50.0 


72 




27.0 


0.60 


This study 




Ireland, mothers 


7 




18.9 


19 




51.4 


11 




29.7 


Ref. 


O'Leary et al., 2006 




Ireland, mothers 


65 




25.0 


119 




45.0 


80 




30.0 


0.70 


This study 




Ireland, fathers 


7 




24.1 


16 




55.2 


6 




20.7 


Ref. 


O'Leary et al., 2006 




Ireland, fathers 


58 




22.0 


127 




48.0 


79 




30.0 


0.58 



simple genetic model, we would expect that relatives would share 
the same risk genotype as the affected individuals and/or moth- 
ers. However, we did not observe this, but instead found that 
maternal relatives had an over-abundance of any risk genotype, 
rather than one specific genotype. We found that more mater- 
nal relatives had one risk genotype, while more paternal relatives 
had no risk genotypes. However, if this were the main explana- 
tion behind the increased number of affected individuals on the 
maternal side, we would expect that relatives more closely related 
to the affected individuals would share more genetic risk factors, 



and thus the prevalence of higher risk genotypes would be greater 
in first-degree relatives than in more distant ones. However, with 
the exception of MTHFR 1298, we did not observe consistent pat- 
terns of increasing numbers of risk genotypes with closer degree 
of relationship to the affected individual. 

The mother's genotype during pregnancy plays a significant 
role in NTD risk, best studied in MTHFR (Yan et al, 2012), and 
the inheritance of these alleles would necessarily come from the 
maternal line. In this case, mothers and possibly affected indi- 
viduals might be expected to have an increased prevalence of 
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Table 4 | Distribution of genotypes and alleles according to relationship within NTD-affected families and lineage, northeast region of the 
Republic of Ireland and Northern Ireland, 1995-2009. 



DHFR Genotype frequencies P b Allele frequencies P b 





Del/Del 


Del/WT 


WT/WT 






Del 






WT 








N 


% 


N 


% 


N 


% 




N 




% 


N 




% 




Affected individual 


z 


I 0.4 


D 


AO. 1 
40.Z 


O 


QQ A 


0.74 


10 




38.5 


16 




61.5 


0.45 


Mother 


7 


18.9 


22 


59.5 


8 


21.6 


0.66 


36 




48.6 


38 




51.4 


0.58 


Father 


6 


20.7 


15 


51.7 


8 


27.6 


0.97 


27 




46.6 


31 




53.4 


0.89 


Sibling 


I O 


ZZ.U 


oO 


oy.o 


I I 


1Q 7 
I O. / 


0.23 


61 




51.7 


57 




48.3 


0.13 


Uncle/aunt 


16 


24.2 


30 


45.5 


20 


30.3 


0.30 


62 




47.0 


70 




53.0 


0.71 


First cousin 


13 


12.4 


61 


58.1 


31 


29.5 


0.11 


87 




41.4 


123 




58.6 


0.13 


Other rel. 


4 


30.8 


3 


23.0 


6 


46.2 


0.06 


11 




42.3 


15 




57.7 


0.73 


1st degree relatives 


26 


20.5 


72 


56.7 


29 


22.8 


0.36 


124 




48.8 


130 




51.1 


0.23 


2nd degree relatives 


17 


24.3 


33 


47.1 


20 


28.6 


0.36 


67 




47.9 


73 




52.1 


0.61 


3rd degree or higher relatives 


16 


14.3 


61 


54.5 


35 


31.2 


0.20 


93 




41.5 


131 




58.5 


0.10 


Paternal line 3 


9 


13.8 


37 


56.9 


19 


29.3 


0.40 


55 




42.3 


75 




57.7 


0.65 


Maternal line 3 


24 


20.9 


55 


47.8 


36 


31.3 




103 




44.8 


127 




55.2 




MTHFD1 




G/G 




G/A 




A/A 






G 






A 






Affected individual 


2 


15.4 


8 


61.5 


3 


23.1 


0.24 


12 




46.2 


14 




53.8 


0.14 


Mother 


12 


32.4 


17 


46.0 


8 


21.6 


0.59 


41 




55.5 


33 




44.6 


0.40 


Father 


14 


48.3 


10 


34.5 


5 


17.2 


0.26 


38 




65.5 


20 




34.5 


0.36 


Sibling 


I D 


z/.z 


oz 


o4.z 


I I 


1Q C 

1 O.D 


0.31 


64 




54.2 


54 




45.8 


0.16 


Uncle/aunt 


34 


51.5 


28 


42.4 


4 


6.1 


0.003 


96 




72.7 


36 




27.3 


0.0008 


First cousin 


34 


32.4 


53 


50.5 


18 


17.1 


0.68 


121 




57.6 


89 




42.4 


0.40 


Other rel. 


3 


23.1 


8 


61.5 


2 


15.4 


0.55 


14 




53.8 


12 




46.2 


0.52 


1st degree relatives 


43 


33.9 


60 


47.2 


24 


18.9 


0.36 


146 




57.5 


108 




42.5 


0.20 


2nd degree relatives 


35 


50.0 


30 


42.9 


5 


7.1 


0.01 


100 




71.4 


40 




28.6 


0.003 


3rd degree or higher relatives 


35 


31.2 


58 


51.8 


19 


17.0 


0.34 


128 




55.4 


96 




42.9 


0.20 


Paternal line 


23 


35.4 


36 


55.4 


6 


9.2 


0.24 


82 




63.1 


48 




36.9 


0.93 


Maternal line 


47 


40.9 


50 


43.5 


18 


15.6 




144 




62.6 


86 




37.4 




MTHFR 1298 




A/A 




A/C 




C/C 






A 






C 






Affected individual 


o 


AR O 

4o.z 


O 


AO. O 
4o.Z 


1 


/.o 


1.00 


18 




69.2 


8 




30.8 


0.88 


Mother 


15 


40.5 


17 


46.0 


5 


13.5 


0.61 


47 




63.5 


27 




36.5 


0.40 


Father 


11 


37.9 


15 


51.8 


3 


10.3 


0.72 


37 




63.8 


21 




36.2 


0.49 


Sibling 


1 Q 
I O 


Oft C 


Jo 




0 


1 u.z 


0.04 


71 




60.2 


47 




39.8 


0.05 


Uncle/aunt 


33 


50.0 


26 


39.4 


7 


10.6 


0.52 


92 




69.7 


40 




30.3 


0.61 


First cousin 


54 


51.4 


44 


41.9 


7 


6.7 


0.22 


152 




72.4 


58 




27.6 


0.09 


Other rel. 


8 


61.5 


4 


30.8 


1 


7.7 


0.49 


20 




76.9 


6 




23.1 


0.31 


1st degree relatives 


44 


34.6 


68 


53.6 


15 


11.8 


0.009 


156 




61.4 


98 




38.6 


0.005 


2nd degree relatives 


35 


50.0 


28 


40.0 


7 


10.0 


0.56 


98 




70.0 


42 




30.0 


0.53 


3rd degree or higher relatives 


60 


53.6 


45 


40.2 


7 


6.2 


0.05 


165 




73.7 


59 




26.3 


0.02 


Paternal line 


30 


46.3 


29 


44.5 


6 


9.2 


0.53 


89 




68.5 


41 




31.5 


0.27 


Maternal line 


63 


54.8 


44 


38.2 


8 


7.0 




170 




73.9 


60 




26.1 




MTHFR 677 




C/C 




C/T 




T/T 






C 






T 






Affected individual 


2 


15.4 


9 


69.2 


2 


15.4 


0.09 


13 




50.0 


13 




50.0 


0.08 


Mother 


17 


46.0 


18 


48.6 


2 


5.4 


0.59 


52 




70.3 


22 




29.7 


0.39 


Father 


11 


37.9 


15 


51.7 


3 


10.4 


0.91 


37 




63.8 


21 




36.2 


0.73 


Sibling 


22 


37.3 


34 


57.6 


3 


5.1 


0.19 


78 




66.1 


40 




33.9 


0.95 


Uncle/aunt 


32 


48.5 


28 


42.4 


6 


9.1 


0.44 


92 




69.7 


40 




30.3 


0.29 


First cousin 


44 


41.9 


47 


44.8 


14 


13.3 


0.33 


135 




64.3 


75 




35.7 


0.56 



(Continued) 
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Table 4 | Continued 
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pb 
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C 


T 








N 


% 


N 


% 


N 


% 




N 


% 


N 


% 




Other rel. 


6 


46.2 


5 


38.5 


2 


15.3 


0.60 


17 


65.4 


9 


34.6 


0.96 


1st degree relatives 


5Z 


a n q 


67 


52.8 


O 
O 


0. 0 


n 1 a 
U. 14 


1 / 1 


67.3 


83 


QO "7 

dz. / 


U. / Z 


2nd degree relatives 


34 


48.6 


30 


42.9 


6 


8.5 


0.53 


98 


70.0 


42 


30.0 


0.32 


3rd degree or higher relatives 


46 


41.1 


50 


44.6 


16 


14.3 


0.12 


142 


63.4 


82 


36.6 


0.22 


Paternal line 


33 


50.8 


26 


40.0 


6 


9.2 


0.38 


92 


70.8 


38 


29.2 


0.16 


Maternal line 


47 


40.9 


52 


45.2 


16 


13.9 




146 


63.5 


84 


36.5 




SLC19A1 




A/A 




A/G 




G/G 






A 


G 






Affected individual 


4 


30.8 


7 


53.9 


2 


15.3 


0.32 


1 5 


57.7 


11 


42.3 


0.15 


Mother 


7 


18.9 


19 


51.4 


1 1 


29.7 


0.99 


33 


44.6 


41 


55.4 


0.93 


Father 


7 


24. 1 


16 


55.2 


6 


20.7 


0.45 


30 


51.7 


28 


48.3 


0.22 


Sibling 


11 


18.6 


35 


59.4 


13 


22.0 


0.34 


57 


48.3 


61 


51.7 


0.31 


Uncle/aunt 


11 


16.7 


27 


40.9 


28 


42.4 


0.04 


49 


37.1 


83 


62.9 


0.07 


First cousin 


16 


15.2 


56 


53.4 


33 


31.4 


0.66 


88 


41.9 


122 


58.1 


0.44 


Other rel. 


2 


15.4 


8 


61.5 


3 


23.1 


0.93 


12 


46.2 


14 


53.8 


0.83 


1st degree relatives 


25 


19.7 


72 


56.7 


30 


23.6 


0.09 


122 


48.0 


132 


52.0 


0.06 


2nd degree relatives 


11 


15.7 


30 


42.9 


29 


41.4 


0.07 


52 


37.1 


88 


62.9 


0.08 


3rd degree or higher relatives 


18 


16.0 


59 


52.7 


35 


31.3 


0.88 


95 


42.4 


129 


57.6 


0.67 


Paternal line 


12 


18.5 


34 


52.3 


19 


29.2 


0.40 


58 


44.6 


72 


55.4 


0.21 


Maternal line 


17 


14.8 


53 


46.1 


45 


39.1 




87 


37.8 


143 


62.2 





£ 'Paternal line consists of uncles and aunts, grandparents, first cousins and first cousins once removed who are related to the affected individual through the father; 
maternal line consists of uncles and aunts, grandparents, first cousins and first cousins once removed who are related to the affected individual through the mother. 
b For specific relatives, each category was compared with all other relatives to calculate P-values. For example, fathers were compared with all other relatives to 
calculate the P-value. For degree of relationship to the affected individual, each category was compared with a combination of the other two to calculate the P-value. 
For example, 1st degree relatives were compared with 2nd and 3rd degree relatives combined to calculate the P-value. P-values were calculated using a chi-square 
test unless one or more cells contained less than 5, in which case Fisher's exact test was used. 



Table 5 | Number of risk genotypes by relationship to affected individuals with NTDs and lineage, northeast region of the Republic of Ireland 
and Northern Ireland, 1995-2009. 



Number of risk genotypes/low metabolizer genotypes P-value a (Fisher's exact test) 
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0 

% 


N 


1 

% 


N 


2 

% 


N 


3 

% 




Maternal line 


24 


20.9 


63 


54.8 


24 


20.9 


4 


3.4 


0.02 


Paternal line 


27 


41.5 


23 


35.4 


12 


18.5 


3 


4.6 




Affected individual or mother 


22 


44.0 


11 


22.0 


15 


30.0 


2 


4.0 


0.004 


All other relatives 


86 


31.6 


130 


47.8 


47 


17.3 


9 


3.3 




Mother 


17 


46.0 


7 


18.9 


12 


32.4 


1 


2.7 


0.005 


All other relatives 


91 


31.9 


134 


47.1 


50 


17.5 


10 


3.5 





3 P-values compare maternal vs. paternal line, affected individual or mother vs. all other relatives, and mother vs. all other relatives. 



genotypes associated with greater NTD risk. In our study, fewer 
affected individuals had the higher folate-metabolizing MTHFR 
677 CC genotype, consistent with earlier reports that the homozy- 
gous TT genotype is a risk factor for NTDs (Shields et al., 
1999; Botto and Yang, 2000). Likewise, affected individuals and 
mothers, as well as mothers only, were more likely to have 2 



risk genotypes. Our study of these associations was limited by 
our small sample size and was not the main objective of this 
study. While our findings were consistent with previous litera- 
ture, these associations alone were not strong enough to explain 
the increased number of affected relatives on the mother's side, 
suggesting that other mechanisms are also at work. 
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Whenever there is preferential inheritance of a condition from 
the maternal vs. paternal line, epigenetic effects are strongly 
suspected. The genes included in our study are not located in 
known imprinted regions. However, folate metabolism is critical 
to the production of one-carbon subunits needed for epigenetic 
modification. Certain variants of folate metabolism genes (pos- 
sibly including those in our study), which are associated with 
NTD risk, might affect proper imprinting and establishment 
or maintenance of other epigenetic changes. These epigenetic 
effects could occur in the maternal genome and would have 
been initiated when the mother was a fetus herself, but may 
also have taken place in the genomes of the mother's parents 
(Crider et al, 2012). If this mechanism were at work, the maternal 
and paternal grandparents' genotypes would be of importance. 
Unfortunately, we were unable to obtain samples from many 
grandparents. However, aunts, uncles, and first cousins provide 
an indirect measure of the grandparents' genotypes, and sig- 
nificant associations for these relatives were seen for DHFR, 
MTHFD1, and SLC19A1. The increased number of risk genotypes 
in maternal vs. paternal relatives is consistent with an epigenetic 
mechanism. 

Our study is the first to look at the genotypes of second- and 
third-degree or more relatives in families affected by NTDs and to 
begin to explore the genetic and epigenetic mechanisms behind 
the increased NTD and overall birth defects risks in these fam- 
ilies. However, our study is preliminary and has a number of 
limitations. The limitations of our study include a relatively small 
sample size, lack of triads (both parents and offspring) among 
distant relatives to allow for heritability testing, and inadequate 
number of samples from grandparents to evaluate epigenetic 
effects. Additionally, more maternal than paternal relatives par- 
ticipated in the study. However, while this might have affected the 
stability of the point estimate for the association of risk genotypes 
with the maternal vs. paternal line, it did not affect the magni- 
tude or direction. The response rate was relatively low for our 
study, and the participation of families and of relatives within 
families could be biased and not reflect the true underlying dis- 
tribution of families with NTDs in Ireland. However, selection 
for and participation in the study was independent of genotype 
and thus participants can be considered close to a random sample 
of relatives. The candidate gene approach used here evaluated a 
small number of polymorphisms and would have missed associ- 
ations with other variants in folate metabolism pathway genes. 
An important limitation of our study (and others) is that we 
were not able to measure the mother's intake of folate and folic 
acid during pregnancy. The genes included in this study function 
in folate metabolism, so that their effect depends on the levels 
of folate present in the body. For example, high intake of folic 
acid and folate could mask genetic effects that would be seen at 
lower intake levels. Lastly, the number of multiple comparisons 
evaluated in our study likely produced some spurious findings. 

Our results are consistent with both genetic and epigenetic 
mechanisms contributing to the increased NTD and overall birth 
defects risk in the maternal lineage. The diversity of genetic risk 
seen in these families suggests that etiology of NTDs is multi- 
faceted. Future studies will benefit from including other relatives, 
especially grandparents, in addition to mothers, fathers, and 



affected individuals, especially if epigenetic mechanisms are to be 
addressed. 

ACKNOWLEDGMENTS 

We are grateful to the families who chose to participate in our 
studies, to Rebecca Lawler, Sharon McGinty, and the volunteers 
who helped to conduct this study, to Claymon Laboratories for 
biochemical tests, to Stanimila Nikolova for technical assistance 
and to Mary E. Cogswell, Quanhe Yang, and Ramal Moonesinghe 
for their insightful comments. The Joseph E. and Marjorie B. 
Jones Foundation of Washington, DC provided partial support to 
Julianne Byrne. No additional external funding was received for 
this study. 

DISCLAIMER 

The findings and conclusions in this report are those of the 
authors and do not necessarily represent the official position of 
the Centers for Disease Control and Prevention. 

REFERENCES 

Botto, L. D., Lisi, A., Bower, C, Canfield, M. A., Dattani, N., De Vigan, C„ et al. 
(2006). Trends of selected malformations in relation to folic acid recommenda- 
tions and fortification: an international assessment. Birth Defects Res. A. Clin. 
Mol. Teratol. 76, 693-705. doi: 10.1002/bdra.20307 

Botto, L. D., and Yang, Q. (2000). 5,10-Methylenetetrahydrofolate reductase gene 
variants and congenital anomalies: a HuGE review. Am. J. Epidemiol. 151, 
862-877. doi: 10.1093/oxfordjournals.aje.a010290 

Brody, L. C, Conley, M., Cox, C, Kirke, P. N., McKeever, M. P., Mills, J. 
L., et al. (2002). A polymorphism, R653Q, in the trifunctional enzyme 
methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohy- 
drolase/formyltetrahydrofolate synthetase is a maternal genetic risk factor for 
neural tube defects: report of the Birth Defects Research Group. Am. J. Hum. 
Genet. 71, 1207-1215. doi: 10.1086/344213 

Byrne, J. (2008). Birth defects in uncles and aunts from Irish families with 
neural tube defects. Birth Defects Res. A. Clin. Mol. Teratol. 82, 8-15. doi: 
10.1002/bdra.20406 

Byrne, J. (2010). Birth defects among maternal first cousins in Irish families with a 
neural tube defect. Ir. J. Med. Sci. 179, 375-380. doi: 10.1007/sll845-009-0381-x 

Byrne, J. (2011). Three generations of matrilineal excess of birth defects in 
Irish families with neural tube defects. Ir. J. Med. Sci. 180, 69-72. doi: 
10.1007/sll845-010-0632-x 

Byrne, J., Cama, A., Reilly, M., Vigliarolo, M., Levato, L., Boni, L., et al. 
(1996). Multigeneration maternal transmission in Italian families with neu- 
ral tube defects. Am. J. Med. Genet. 66, 303-310. doi: 10.1002/(SICI)1096- 
8628(19961218)66:3<303::AID-AJMG13>3.3.CO;2-2 

Carter, T. C, Pangilinan, P., Troendle, J. F., Molloy, A. M., VanderMeer, J., Mitchell, 
A., et al. (2011). Evaluation of 64 candidate single nucleotide polymorphisms as 
risk factors for neural tube defects in a large Irish study population. Am. J. Med. 
Genet. A 155, 14-21. doi: 10.1002/ajmg.a.33755 

Chango, A., Emery-Fillon, N., de Courcy, G. P., Lambert, D., Pfister, M., Rosenblatt, 
D. S., et al. (2000). A polymorphism (80G->A) in the reduced folate carrier gene 
and its associations with folate status and homocysteinemia. Mol. Genet. Metab. 
70, 310-315. doi: 10.1006/mgme.2000.3034 

Christensen, K. E., Rohlicek, C. V., Andelfinger, G. U., Michaud, J., Bigras, J. L., 
Richter, A., et al. (2009). The MTHFD1 p.Arg653Gln variant alters enzyme 
function and increases risk for congenital heart defects. Hum. Mutat. 30, 
212-220. doi: 10.1002/humu.20830 

Coffey, V. P. (1983). Neural tube defects in Dublin 1953-1954 and 1961-1982. Ir. 
Med. J. 76,411-413. 

Crider, K. S., Yang, T. P., Berry, R. J., and Bailey, L. B. (2012). Folate and DNA 
methylation: a review of molecular mechanisms and the evidence for folate's 
role. Adv. Nutr. 3, 21-38. doi: 10.3945/an.lll.000992 

Czeizel, A. E., and Dudas, I. (1992). Prevention of the first occurrence of neural- 
tube defects by periconceptional vitamin supplementation. N. Engl. J. Med. 327, 
1832-1835. doi: 10.1056/NEJM199212243272602 



Frontiers in Genetics | Epigenomics and Epigenetics 



November 2013 | Volume 4 | Article 223 | 8 



Fisk Green et al. 



Folate gene variants in families 



Deak, K. L., Siegel, D. G., George, T. M., Gregory, S., Ashley-Koch, A., Speer, M. C, 
et al. (2008). Further evidence for a maternal genetic effect and a sex-influenced 
effect contributing to risk for human neural tube defects. Birth Defects Res. A. 
Clin. Mol. Teratol 82, 662-669. doi: 10.1002/bdra.20511 

De Marco, P., Calevo, M. G., Moroni, A., Merello, E., Raso, A., Finnell, R. H., 
et al. (2003). Reduced folate carrier polymorphism (80A->G) and neural tube 
defects. Eur. J. Hum. Genet. 11, 245-252. doi: 10.1038/sj.ejhg.5200946 

Frosst, P., Blom, H. J., Milos, R., Goyette, P., Sheppard, C. A, Matthews, R. G., 
et al. (1995). A candidate genetic risk factor for vascular disease: a common 
mutation in methylenetetrahydrofolate reductase. Nat. Genet. 10, 1 1 1-1 13. doi: 
10.1038/ng0595-lll 

Kalmbach, R. D., Choumenkovitch, S. F., Troen, A. P., Jacques, P. F., D'Agostino, 
R., and Selhub, J. (2008). A 19-base pair deletion polymorphism in dihy- 
drofolate reductase is associated with increased unmetabolized folic acid in 
plasma and decreased red blood cell folate. /. Nutr. 138, 2323-2327. doi: 
10.3945/jn.l08.096404 

Kirke, P. N., Mills, J. L., Molloy, A. M., Brody, L. C, O'Leary, V. B., 
Daly, L., et al. (2004). Impact of the MTHFR C677T polymorphism on 
risk of neural tube defects: case-control study. BMJ 328, 1535-1536. doi: 
10.1136/bmj.38036.646030.EE 

Mills, J. L., Molloy, A. M., Parle-McDermott, A., Troendle, J. E, Brody, L. C, Conley, 
M. R., et al. (2008). Folate-related gene polymorphisms as risk factors for cleft 
lip and cleft palate. Birth Defects Res. A. Clin. Mol. Teratol. 82, 636-643. doi: 
10.1002/bdra.20491 

Molloy, A. M., Brody, L. C, Mills, J. L., Scott, J. M., and Kirke, P. N. (2009). The 
search for genetic polymorphisms in the homocysteine/folate pathway that con- 
tribute to the etiology of human neural tube defects. Birth Defects Res. A. Clin. 
Mol. Teratol. 85, 285-294. doi: 10.1002/bdra.20566 

MRC Vitamin Study Research Group. (1991). Prevention of neural tube defects: 
results of the Medical Research Council Vitamin Study. Lancet 338, 131-137. 
doi: 10.1016/0140-6736(91)90133-A 

Mynett-Johnson, L. A., Keenan, C, Black, I. L., Livingstone, W. J., Lawler, M., 
Roche, H. M., et al. (2002). Thermolabile methylenetetrahydrofolate reductase 
(C677T): frequency in the Irish population. Ir. J. Med. Sci. 171, 37-39. doi: 
10.1007/BF03168940 

O'Connell, J. R., and Weeks, D. E. (1998). PedCheck: a program for identification of 
genotype incompatibilities in linkage analysis. Am. ]. Hum. Genet. 63, 259-266. 
doi: 10.1086/301904 

O'Leary, V. B., Pangilinan, E, Cox, C, Parle-McDermott, A., Conley, M., Molloy, A. 
M., et al. (2006). Reduced folate carrier polymorphisms and neural tube defect 
risk. Mol. Genet. Metab. 87, 364-369. doi: 10.1016/j.ymgme.2005.09.024 

Pangilinan, E, Molloy, A. M., Mills, J. L., Troendle, J. F., Parle-McDermott, A., 
Signore, C, et al. (2012). Evaluation of common genetic variants in 82 can- 
didate genes as risk factors for neural tube defects. BMC Med. Genet. 13:62. doi: 
10.1186/1471-2350-13-62 

Parle-McDermott, A., Kirke, P. N., Mills, J. L., Molloy, A. M., Cox, C, O'Leary, 
V. B., et al. (2006). Confirmation of the R653Q polymorphism of the tri- 
functional Cl-synthase enzyme as a maternal risk for neural tube defects in 
the Irish population. Eur. J. Hum. Genet. 14, 768-772. doi: 10.1038/sj.ejhg. 
5201603 

Parle-McDermott, A., Mills, J. L., Kirke, P. N., O'Leary, V. B., Swanson, D. A., 
Pangilinan, E, et al. (2003). Analysis of the MTHFR 1298A->C and 677C->T 
polymorphisms as risk factors for neural tube defects. /. Hum. Genet. 48, 
190-193. doi: 10.1007/sl0038-003-0008-4 

Parle-McDermott, A., Pangilinan, E, Mills, J. L., Kirke, P. N., Gibney, E. R., 
Troendle, J., et al. (2007). The 19-bp deletion polymorphism in intron-1 of 



dihydrofolate reductase (DHFR) may decrease rather than increase risk for 
spina bifida in the Irish population. Am. J. Med. Genet. A 143A, 1174-1180. 
doi: 10.1002/ajmg.a.31725 

Shaw, G. M., Lu, W., Zhu, H., Yang, W., Briggs, F. B., Carmichael, S. L., et al. (2009). 
118 SNPs of folate-related genes and risks of spina bifida and conotruncal heart 
defects. BMC Med. Genet. 10:49. doi: 10.1186/1471-2350-10-49 

Shields, D. C, Kirke, P. N., Mills, J. L., Ramsbottom, D., Molloy, A. M., Burke, EL, 
et al. (1999). The "thermolabile" variant of methylenetetrahydrofolate reductase 
and neural tube defects: an evaluation of genetic risk and the relative impor- 
tance of the genotypes of the embryo and the mother. Am. }. Hum. Genet. 64, 
1045-1055. doi: 10.1086/302310 

Stanislawska-Sachadyn, A., Brown, K. S., Mitchell, L. E., Woodside, J. V, Young, 
I. S., Scott, J. M., et al. (2008). An insertion/deletion polymorphism of 
the dihydrofolate reductase (DHFR) gene is associated with serum and red 
blood cell folate concentrations in women. Hum. Genet. 123, 289-295. doi: 
10.1007/s00439-008-0475-y 

Sutton, M., Daly, L. E., and Kirke, P. N. (2008). Survival and disability in a cohort 
of neural tube defect births in Dublin, Ireland. Birth Defects Res. A. Clin. Mol. 
Teratol. 82, 701-709. doi: 10.1002/bdra.20498 

van der Put, N. M., Gabreels, F., Stevens, E. M., Smeitink, J. A., Trijbels, F. J., Eskes, 
T. K., et al. (1998). A second common mutation in the methylenetetrahydro- 
folate reductase gene: an additional risk factor for neural-tube defects? Am. J. 
Hum. Genet. 62, 1044-1051. doi: 10.1086/301825 

Wang, X. W., Luo, Y. L., Wang, W, Zhang, Y, Chen, Q., and Cheng, Y. L. 
(2012a). Association between MTHFR A1298C polymorphism and neural tube 
defect susceptibility: a metaanalysis. Am. J. Ohstet. Gynecol 206, 251.el-7. doi: 
10.1016/j.ajog.2011. 12.021 

Wang, H. G, Wang, J. L., Zhang, J., Zhao, L. X., Zhai, G. X., Xiang, Y. 
Z., et al. (2012b). Reduced folate carrier A80G polymorphism and sus- 
ceptibility to neural tube defects: a meta-analysis. Gene 510, 180-184. doi: 
10.1016/j.gene.2012.02.020 

Whitehead, A. S., Gallagher, P., Mills, J. L., Kirke, P. N., Burke, H., Molloy, A. M., 
et al. (1995). A genetic defect in 5,10 methylenetetrahydrofolate reductase in 
neural tube defects. QJM 88, 763-766. 

Yan, L., Zhao, L., Long, Y, Zou, P., Ji, G., Gu, A., et al. (2012). Association of 
the maternal MTHFR C677T polymorphism with susceptibility to neural 
tube defects in offsprings: evidence from 25 case-control studies. PLoS ONE 
7:e41689. doi: 10.1371/journal.pone.0041689 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 09 July 2013; accepted: 15 October 2013; published online: 06 November 
2013. 

Citation: Fisk Green R, Byrne J, Crider KS, Gallagher M, Koontz D and Berry R] 
(2013) Folate-related gene variants in Irish families affected by neural tube defects. 
Front. Genet. 4:223. doi: 10.3389/fgene.2013.00223 

This article was submitted to Epigenomics and Epigenetics, a section of the journal 
Frontiers in Genetics. 

Copyright © 2013 Fisk Green, Byrne, Crider, Gallagher, Koontz and Berry. This is an 
open-access article distributed under the terms of the Creative Commons Attribution 
License (CC BY). The use, distribution or reproduction in other forums is permitted, 
provided the original author(s) or licensor are credited and that the original publica- 
tion in this journal is cited, in accordance with accepted academic practice. No use, 
distribution or reproduction is permitted which does not comply with these terms. 



www.frontiersin.org 



November 2013 | Volume 4 | Article 223 | 9 



